arXiv:1504.03204vl [astro-ph.SR] 13 Apr 2015 


Astronomy & Astrophysics manuscript no. AGCar_preprint ©ESO 2015 

April 14, 2015 


The Herschel view of the nebula around the luminous 
blue variable star AG Carinae 

C. Vamvatira-Nakou* ** , D. Hutsemekers'■***, P. Royer^, N.L.J. Cox^, Y. Naze*’****, G. Rauw’, 

C. Waelkens^, and M. A. T. Groenewegen^ 

' Institut d’Astrophysique et de Geophysique, Universite de Liege, Allee du 6 aout, 17 - Bat. B5c, B-4000 Liege, Belgium 
e-mail: vamvatiraSastro. ulg .ac.be 

^ Instituut voor Sterrenkunde, KU Leuven, Celestijnenlaan 200D, Bus 2401, B-3001 Leuven, Belgium 
^ Koninklijke Sterrenwacht van Belgie, Ringlaan 3, B-1180 Brussels, Belgium 

Preprint online version: April 13, 2015 


ABSTRACT 

Far-infrared Herschel PACS imaging and spectroscopic observations of the nebula around the luminous blue variable (LBV) star 
AG Car have been obtained along with optical imaging in the Hcr-r [N ii] filter. In the infrared light, the nebula appears as a clumpy 
ring shell that extends up to 1.2 pc with an inner radius of 0.4 pc. It coincides with the Ho- nebula, but extends further out. Dust 
modeling of the nebula was performed and indicates the presence of large grains. The dust mass is estimated to be ~ 0.2 Mo. The 
infrared spectrum of the nebula consists of forbidden emission lines over a dust continuum. Apart from ionized gas, these lines also 
indicate the existence of neutral gas in a photodissociation region that surrounds the ionized region. The abundance ratios point 
towards enrichment by processed material. The total mass of the nebula ejected from the central star amounts to ~ 15 Mq, assuming a 
dust-to-gas ratio typical of LB Vs. The abundances and the mass-loss rate were used to constrain the evolutionary path of the central 
star and the epoch at which the nebula was ejected, with the help of available evolutionary models. This suggests an ejection during a 
cool LBV phase for a star of ~ 55 Mq with little rotation. 
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1. Introduction 


The term “luminous blue variable” (LBV) was used by Conti 
111984b for the first time and referred to hot luminous massive 
variable stars that are evolved, but are not Wolf-Rayet (WR) 
stars. Nowadays, LB Vs, also known as S Doradus variables, 
are considered to be massive evolved stars mainly character¬ 
ized by a) high luminosity, ~ 10® Lq; b) photometric vari¬ 
ability with amplitudes from ~0.1 mag (small oscillations) up 
to > 2 mag (giant eruptions); and c) high mass-loss rates ~ 
10^® - 10 Mq yr ' (Humphreys & Davidson ll994b . Their lo¬ 
cation on the Hertzsprung- Russell (HR) diagram is in the upper 
left part although some of them undergo occasional excursions 
to the right. 

An early-type O star with initial mass > 30 Mq evolves to a 
WR star by losing a significant fraction of its initial mass. Lu¬ 
minous blue variables represent a short stage (~ 10^ - 10® yr) 
in this evolutionary path according to current evolutionary sce¬ 
narios (Maeder & Mevnet l2010b . Although stellar winds can be 
responsible for stellar mass-loss, the mass-loss rates of O stars 
have been revised downward in the past few years (Bouret et 
al. 120051 Lullerton et al. 120061 Puls et al. 12008b . Consequently, 
episodes of extreme mass-loss during an intermediate evolution¬ 
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ary phase, like a LBV or a red supergiant (RSG) phase, are now 
thought to play a key role, which is why the study of the LB Vs 
and their circumstellar environments is crucial for understanding 
massive star evolution. 

Such extreme mass-loss leads to the formation of ejected 
nebulae, which have been observed around many LB Vs (Hut- 
semekers ll994l Nota et al. 11995b . They are classified into three 
categories according to their morphology: shell nebulae, fila¬ 
mentary nebulae and peculiar morphologies (Nota et al. 11995b . 
The study of these circumstellar environments can reveal the 
mass-loss history of the central star since they are formed by 
the material that has been ejected from the central star in a previ¬ 
ous evolutionary phase. Dust and molecular gas (CO) have been 
revealed by infrared and millimeter studies of LBV nebulae (Mc¬ 
Gregor et al. 119881 Hutsemekers 1 1 9971 Nota et al. 12002b . 

Some LB Vs are surrounded by more than one nebulae. This 
is the case of the LBV G79.29 h- 0.46. Near-infrared and mil¬ 
limeter data analyzed by Jimenez-Esteban et al. (12010b revealed 
multiple shells around this star. Infrared observations by the 
Herschel Space Observatory (Pilbratt et al. 12010b revealed a sec¬ 
ond nebula around the LBV WRAY 15-751 (Vamvatira-Nakou 
et al. l20T3l l 

AG Car (=HD 94910 MRAS 10541-6011) is a well-studied 
prototypical LBV. Its variability was first discovered by Wood 
( 11914b . It was first classified as a P Cygni star by Cannon ( 11916b 
and finally classified as a LBV by Humphreys (11989b . Numer¬ 
ous studies show that this star exhibits photometric and spectro¬ 
scopic variability. In the optical V-band, the photometric changes 
during the S Dor cycle are about 2 mag on a timescale of 5- 
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10 years (Stahl [198^ van Genderen et al. 119881 Leitherer et al. 
119921 Sterken et al. 119961) . In addition, smaller variations of 0.1- 
0.5 mag on a timescale of about 1 year were discovered (van 
Genderen et al. 11997b . During the periods of visual minimum, 
AG Car has a spectrum of a WR star, with Ofpe/WN9 spectral 
type according to Stahl (11986b and with WN11 according to a 
more recent study by Smith et al. (11994b . During the periods of 
visual maximum, AG Car’s spectrum corresponds to that of an 
early-A hypergiant (Wolf & Stahl 1 198^ Stahl et al. 12001b . 

Humphreys et al. (11989b concluded that the distance to AG 
Car is 6 + 1 kpc, based on the calculated kinematic distance 
and on the observed variation of the interstellar extinction with 
distance. This result was confirmed by Hoekzema et al. (11992b . 
based again on the extinction versus distance relation. Stahl et 
al. (12001b suggested a slightly lower distance of 5-6 kpc based 
on their calculations of the heliocentric systemic velocity of 
AG Car (10 + 5 km s“'), which is compatible with the value 
of Humphreys et al. (11989b considering the errors. Groh et al. 
(12009b calculated a similar systemic velocity. Consequently, the 
value of 6 + 1 kpc that encompasses all measurements is adopted 
for all calculations in this study. 

Lamers et al. (11989b calculated logL/L© = 6.2 + 0.2 for the 
luminosity of AG Car and showed that it remains constant during 
the light variations of the star, as was also found for other LB Vs 
(R71: Wolf et a l.flMn R127: Stahl and Wolf fT986b . Later on, 
Leitherer et al. (11994b . in their study of the stellar wind of AG 
Car, found a slightly lower bolometric luminosity of 6.0 + 0.2 
based on ultraviolet observations combined with visual and near- 
infrared photometry. They also confirmed the nonvariability of 
the bolometric luminosity during the S Dor cycle. Given such a 
high value of bolometric luminosity, AG Car is well above the 
Humphreys-Davidson limit (Humphreys and Davidson 11979b . 
the limit above which a massive star becomes unstable and high 
mass loss episodes take place. However, in a recent study of the 
fundamental parameters of AG Car during visual minima, Groh 
et al. (12009b concluded that the bolometric luminosity of AG 
Car does change during the S Dor cycle. They obtained a maxi¬ 
mum value of the bolometric luminosity during minimum phase 
of logL/Lo = 6.18, with a variation amplitude of A(logL/Lo) ~ 
0.17 dex. This luminosity variation lies inside the limits of the 
previously calculated values considering the errors. In all these 
studies, the distance 6+1 kpc was used. 

Thackeray (11950b discovered that AG Car is surrounded by 
a nebulous shell that has the shape of an elliptical ring with 
nonuniform intensity. He measured the size of the nebula to be 
39" X 30" and the width of the ring to be about 5". Stahl (11987b 
studied direct CCD imaging data of AG Car and suggested that 
very likely the nebula around this star is the result of a heavy 
mass-loss episode that took place during an S Dor outburst and 
not the result of interstellar material that was swept up by the 
stellar wind. Smith (11991b studied the dynamics of AG Car neb¬ 
ula based on spectroscopic observations in the optical waveband. 
She measured an average nebular expansion velocity of 70 km 
s ' and concluded that the shell expansion is roughly symmetri¬ 
cal. She also reported the presence of a jet-like bipolar mass out¬ 
flow, expanding with a velocity of 83 km s ' and distorting the 
northeastern edge of the shell. Nota et al. (11992b studied the neb¬ 
ula around AG Car using high-resolution coronographic imaging 
and spectroscopic data so as to constrain its geometry. They con¬ 
firmed the value of the expansion velocity of Smith (11991b and 
concluded that this nebular shell shows a deviation from spher¬ 
ical symmetry, based on the observed radial velocity variations 
and on the gas distribution observed in the images. 


Voors et al. (12000b studied AG Car in the infrared waveband 
by modeling ground-based infrared images taken at about 10 pm 
and ISO spectroscopic observations, from which they derived the 
properties of the circumstellar dust. The dust shell is detached 
and slightly elongated. The ionized gas appears co-spatial with 
the dust. Polycyclic aromatic hydrocarbons (PAHs) are present. 
The dust shell contains mostly large grains, although very large 
grains are present and also a population of small, warm, out 
of thermal equilibrium grains that produce continuum and PAH 
bands emission. 

Duncan and White (12002b observed the AG Car nebula at the 
radio wavelengths (3 and 6 cm). Their 3 cm wavelength radio 
image revealed a nebula with a detached ring shape, very similar 
to the morphology in the Ha;+[N ii] filter (see Sect. (3). 

Nota et al. (12002b detected '^CG 7=1^0 and 7 = 2 —» 1 
emission from AG Car for the first time. The CO line profiles in¬ 
dicate a region of molecular gas that is close to the star, expand¬ 
ing slowly and not originating from the gaseous nebula. They 
argued that the most plausible scenario to explain the observed 
profile is the presence of a circumstellar disk. 

Weis (12008b . using deep Her imaging, reported the presence 
of diffuse emission in the form of a cone-like extension to the 
north of the AG Car nebula and concluded that it is clearly part 
of the nebula. It extends up to 28" outside the nebula, increasing 
its size by about two times. It has a higher radial velocity than 
the ring. 

In this study we present an analysis and discussion of the 
images and the spectrum of the AG Car nebula taken by PACS 
(Photodetector Array Camera and Spectrometer, Poglitsch et 
al. 12010b . which is one of the three instruments on board the 
Herschel Space Observatory. The paper is organized as follows. 
In Sect. |3the observations and the data reduction procedure are 
presented. A description of the nebular morphology based on 
these observations is given in Sect. [3 The dust continuum emis¬ 
sion is modeled in Sect. |4] while the emission line spectrum is 
analyzed in Sect. |5] In Sect. |6]a general discussion is presented 
and in Sect. |7]conclusions are drawn. 


2. Observations and data reduction 

2.1. Infrared observations 

The infrared observations include imaging and spectroscopy of 
the AG Car nebula and were carried out using PACS in the 
framework of the Mass-loss of Evolved StarS (MESS) Guaran¬ 
teed Time Key Program (Groenewegen et al. l201 lb . 

The PACS imaging observations were carried out on Au¬ 
gust 12, 2010, which corresponds to HerscheVs observational 
day (OD) 456. The observing mode was the scan map in which 
the telescope slews at constant speed (in our case the “medium” 
speed of 20"/s) along parallel lines in order to cover the re¬ 
quired area of the sky. Two orthogonal scan maps were obtained 
for each filter. Our final data set consists of maps at 70, 100 
and 160 pm. The observation identification numbers (obsID) of 
the four scans are 1342202927, 1342202928, 1342202929, and 
1342202930; each scan has a duration of 157s. 

To perform the data reduction we used the Herschel Inter¬ 
active Processing Environment (HIPE, Ott 12010b up to level 1. 
Subsequently, the data were reduced and combined using the 
Scanamorphos software (Roussel 12013b . The pixel size in the 
final maps is 2" in the blue channel (70, 100 pm) and 3" in the 
red channel (160 pm). It should be mentioned that the Herschel 
PACS point spread function (PSP) full widths at half maximum 
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Fig. 1. Images of the nebula around AG Car. Top: PACS images at 70 yum, 100 //m, and 160 //m from left to right. Bottom: the Htt+[N ii] image 
(left), the continuum image (right) and the image resulting from the subtraction of the continuum image from the Hu+iN ii] image after correcting 
for the position offsets and for the different filter transmissions using field stars (middle). The size of each image is 1.5' x 1.5'. The scale on the 
right corresponds to the surface brightness (arbitrary units). North is up and east is to the left. 



Fig. 2. Left: View of the nebula in the optical. The 
bright Ha+iNii] ring nebula is illustrated in red 
(also shown in Fig. [T] at the same scale) while the 
fainter Ha;+[Nii] emission is shown in blue, dif¬ 
ferently scaled in surface brightness, revealing the 
northern extension (Weis l2008t . The continuum im¬ 
age that outlines the dust scattered nebula is pre¬ 
sented in green. Right: Contour image of the optical 
emission from the nebula (green lines) superposed 
on the infrared image of the nebula at 70 jim (also 
shown in Fig. [T] at the same scale). The size of each 
image is 1.5' x 1.5'. North is up and east is to the 
left. 


(FWHMs) are 5'.'2, 7'.'7, and 12" at 70 fim, 100 fim and 160 fim, 
respectively. 

The spectrum of the AG Car nebula was taken on June 5, 
2010 (OD 387), with the PACS integral-field spectrometer that 
covers the wavelength range from 52 /tm to 220 fim in two 
channels that operate simultaneously in the blue 52-98 fim band 
(second order; B2A 52-73 fim and B2B 70-105 fim), and the 
red 102-220 fim band (first order; RIA 133-220 fim and RIB 
102-203 fim). Simultaneous imaging of a 47" x 47" field of 
view is provided that is resolved in 5 x 5 square spatial pixels 
(i.e., spaxels). The two-dimensional field of view is then rear¬ 
ranged along a 1 X 25 pixel entrance slit for the grating via an 


image sheer employing reflective optics. Its resolving power is 
A/6A ~ 940 - 5500 depending on the wavelength. The observ¬ 
ing template was the spectral energy distribution (SED) that pro¬ 
vides a complete coverage between 52 and 220 fim. The two ob- 
sIDs are 1342197792 and 1342197793. HIPE was also used for 
the data reduction. The standard reduction steps were followed 
and in particular the subtraction of the background spectrum ob¬ 
tained through chopping and nodding. 
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Fig. 3. Three-color image (70 /jm in blue, 100 /jm in green and 160 pm 
in red) of the AG Car nebula. The LBV nebula appears located inside a 
cavity in the interstellar medium. The size of the image is 7' x 7'. North 
is up and east is to the left. 

2.2. Visible observations 

The optical images of the AG Car nebulae were obtained on 
April 6, 1995, with the 3.6 m telescope at the European South¬ 
ern Observatory (ESO), La Silla, Chile. A series of short (Is-10s) 
and longer (30s-60s) exposures were secured in a HaH-[N ii] fil¬ 
ter (Ac - 6560.5A; EWHM =62.2A) and in a continuum filter 
just redwards (Ac = 6644.7A; EWHM =61.0A). The EEOSCl 
camera was used in its coronographic mode: for the longer ex¬ 
posures, the 8" circular coronographic mask was inserted in the 
aperture wheel and positioned on the central star while the Lyot 
stop was inserted in the grism wheel (Melnick et al. ll989l l. The 
frames were bias-corrected and flat-fielded. The CCD pixel size 
was 07605 on the sky. The night was photometric and the seeing 
around 172. In order to properly calibrate the images, three spec- 
trophotometric standard stars and three planetary nebulae with 
known Ha flux were observed in the HaH-[N ii] Alter. 

3. Morphology of the nebula 

The three PACS infrared images of the nebula around the LBV 
AG Car at 70 pm, 100 pm, and 160 pm along with images taken 
at optical wavelengths are shown in Eigs.[T]and |2] 

The infrared images reveal a dusty shell nebula with a 
clumpy ring morphology that is clearly detached from the central 
star (not visible at these wavelengths). The brightness of the neb¬ 
ula is not uniform, with the southwestern part being the brightest 
one. This agrees with the nonuniform brightness of the images at 
12.5 pm and 12.78 pm analyzed by Voors et al. (12000b . An axis 
at position angle (PA) ~ 160° can be defined (north: PA = 0°; 
east: PA = 90°). Cuts of the 70 pm PACS image through PA ~ 
160° and PA ~ 70° (Eig.lHi show peaks at a radius of about 14", 
although with different intensities. The ring extends up to ~ 42". 
These values correspond to 0.4 pc and 1.2 pc, respectively, at a 
distance of 6 kpc. On top of this central circular nebula, there is 


a northern faint extension. No other more extended nebula asso¬ 
ciated with the star can be detected. A three-color infrared image 
of the nebula and its environment is illustrated in Pig. [2 which 
shows best this faint northern extension first reported by Weis 
(12008b . although the resolution is not good enough to reveal its 
detailed structure. The LBV nebula seems to be located in a cav¬ 
ity in the interstellar medium, probably excavated by the star in 
a previous evolutionary phase. The radius of this empty cavity is 
about 2.5', which corresponds to 4.4 pc at a distance of 6 kpc. 

The Han-[Nil] image (Pig. [TJ shows that the gas nebula 
around AG Car forms an elliptical shell with an average outer 
radius of ~ 20" and an inner radius of ~11". To more accu¬ 
rately investigate the morphology of the nebula in the optical, a 
three-color image is presented in Pig. |2] The fainter nebular op¬ 
tical emission reveals the northern extension described in Weis 
(12008b . The nebula that surrounds the bright ring extends up to 
~ 23", while the north extension goes up to ~ 36" from the cen¬ 
ter of the nebula. These numbers correspond to 0.7 pc and 1 pc, 
respectively, at a distance of 6 kpc, in agreement with the size of 
the nebula given in Weis (120081 l2UTTb . 

In the Ha light, the dynamics points to a spherically expand¬ 
ing shell distorted by a more extended bipolar nebula (Smith 
119911 Nota et al. 11992b . In projection on the sky, the shell ap¬ 
pears as an elliptical ring with PA ~ 131°, different from the 
infrared shell PA. Nevertheless, the contour image of the optical 
emission (both bright and faint) superimposed on the infrared 
image of the nebula at 70 pm as illustrated in the right panel of 
Fig.E shows that the overall morphology of the gas nebula is 
similar to the infrared dust morphology, although the Ha-(-[N ii] 
ring nebula appears slightly smaller and more elliptical than the 
infrared nebula. The bright region at the southwestern part of 
the gas nebula coincides with the bright region of the infrared 
dust nebula. The northern faint extended structure, which ap¬ 
pears both in the infrared and in the optical, is likely a lobe of the 
bipolar nebula. The extension seen to the south in the Ha-(-[N ii] 
map and in the velocity maps of Smith (11991b could constitute 
a part of a second fainter lobe. The system may thus consist of 
a typical bipolar nebula seen roughly through the poles (i.e., at 
inclination < 30°), with two faint lobes and a bright waist. 

The nebula is also clearly detected in the optical continuum 
Alter (Pig. [T]i, indicating significant dust scattering. The mor¬ 
phology of the dust reflection nebula appears somewhat different 
from the morphology of the Ha-(-[N ii] emission. This reflection 
nebula around AG Car was first described in Paresce and Nota 
(11989b and its stunning structure resolved with the Hubble Space 
Telescope (Nota et al. 11995111996b . More specifically, in the op¬ 
tical continuum image (Pigs. [T] and El the ring appears circu¬ 
lar and very clumpy with a jet-like feature that starts from the 
central part of the nebula and extends towards the southwestern 
part, which is the brightest region of the dust and gas emission. 
The northern extension of the nebula is not detected in the op¬ 
tical continuum. As does the Ha-i-[Nn] ring nebula, the optical 
continuum emission appears slightly inside the infrared contin¬ 
uum emission, although a detailed comparison is difficult given 
the lower spatial resolution of the PACS images. The difference 
in morphology between the optical continuum and the Ha-(-[N ii] 
nebulae may arise from anisotropic illumination and thus ioniza¬ 
tion of different parts of the nebula owing to its clumpy structure. 

Considering the nebular expansion velocity, Vexp, of 70 km 
s ' measured by Smith (11991b . the kinematic age, fkin, of the 
nebula can be estimated. As mentioned above, the nebula in the 
infrared extends to 1.2 pc in radius, r, so it has a kinematic age of 
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fkin = t-/Vexp = 1.7x10'* yeuTS. The temporal difference between 
the inner and the outer radius of the nebula is 1.1x10'* years. 

4. Dust continuum emission 

Integrated flux densities were derived for the nebular shell at the 
three PACS wavelengths by performing aperture photometry on 
the PACS images. We also used imaging data taken from the 
archives of the Infrared Astronomical Satellite (IRAS) mission 
(Neugebauer et al. ll984t and the Infrared Astronomical Mission 
AKARI (Murakami et al. 120071) . In these archival data we did 
not include the IRAS observation at 12 fim because it is only an 
upper limit or at 100 fim because it is not of high quality. We 
note that the beam size of the IRAS and AKARI observations 
is large enough to fully encompass the ring nebula around AG 
Car. Images from the SPIRE (Spectral and Photometric Imaging 
Receiver, Griffin et al. 120101) instrument on board the Herschel 
Space Observatory were also included. They were taken from the 
observations of the Herschel Infrared Galactic Plane survey (Hi- 
GAL, Molinari et al. 120101) made immediately public for legacy. 
The Hi-GAL observations of the field around AG Car were re¬ 
trieved from the archive processed up to level 2. Only the maps at 
250 jum and at 350 //m were used, because at 500 jum the nebula 
is very faint so that flux determination is highly uncertain. Inte¬ 
grated flux densities were derived for the nebular shell at these 
two SPIRE wavelengths by performing aperture photometry on 
the images. 

We applied photometric color correction to all flux densi¬ 
ties derived from the observations of these three space missions. 
With this correction, the monochromatic flux densities that refer 
to a constant energy spectrum are converted to the true object 
SED flux densities at the photometric reference wavelengths of 
each instrument. We used the flux density ratio to derive the color 
temperature for the color correction of the IRAS data and then 
we chose the corresponding color-correction factor (Beichman 
et al. ll988l) . The ratio R (25,60) corresponds to a temperature of 
140 K, so the correction factors for this temperature were used. 
To color correct the AKARI EIS and IRC data we fitted a black- 
body to the two data sets independently using the 25 jum IRAS 
observation because we needed a measurement near the maxi¬ 
mum of the curve. These fits led us to adopt the color-correction 
factors that correspond to a temperature of 150 K for EIS (Yama- 
mura et al. 120101) and 220 K for IRC data (Lorente et al. 120081 . 
To estimate the color correction of the Herschel-PACS data, we 
fitted a blackbody, considering again the 25 /tm IRAS observa¬ 
tion. This fit gave a temperature of 150 K, for which we adopted 
the corresponding correction factor (Muller et al. 1201II) . Eor the 
SPIRE data, the instructions for color correction given in the 
SPIRE HandbookQ were followed. 

The corrected measurements are presented in Table [1] They 
enabled us to construct the infrared SED of the nebula, along 
with the archived spectrum that was part of the observations 
carried out by the Infrared Space Observatory (ISO) mission 
(Kessler et al. 119961) . A detailed discussion on this ISO-LWS 
spectrum can be found in Voors et al. (120001) . 

The infrared SED of the nebula around AG Car obtained at 
different epochs with the various instruments is shown in Eig.|5] 
All these measurements agree very well within the uncertainties 
except the one at 90 pm. This is likely due to the uncertainty on 
the color correction that was stronger for the 90 pm (AKARI- 
EIS) data point. 


* http://herschel.esac.esa.int/Docs /SPIRE/spire_handbook.pdf 


Table 1. Color-corrected nebular flux densities. 


Spacecraft-Instrument Date 

4 

(pm) 

Fy 

(Jy) 

Error 

(Jy) 

IRAS“ 

1983 

25 

187.5 

9.4 



60 

177.7 

28.4 

AKARI-IRC* 

2007 

9 

9.04 

0.14 



18 

119.2 

0.24 

AKARI-EIS^ 

2007 

65 

229.3 

6.53 



90 

81.0 

3.9 



140 

45.5 

5.2 



160 

35.1 

4.6 

Herschel-PACS"* 

2010 

70 

173 

2 



100 

103 

3 



160 

42 

3 

Herschel-SPIRE" 

2010 

250 

8.1 

2 



350 

3.0 

1 


Notes. Data from: IRAS Point Source Catalog (Beichman et al. 

119881 . Akari/TRC Point Source Catalogue (Ishihara et al. 12010b . 

Akari/FIS Bright Source Catalogue (Yamamura et al. 12010b . This 
work. Observations of Hi-GAL (Molinari et al. 12010b retrieved from 
the Herschel archive. 

A model of the dust nebula around the LBV AG Car has 
been carried out in the past by Voors et al. (12000b . They used a 
one-dimensional radiative transfer code to fit both imaging and 
spectroscopic infrared data. To further constrain the dust proper¬ 
ties, we use the AKARI archive imaging data and the new PACS 
and SPIRE imaging data in addition to the IRAS imaging data 
and the ISO infrared spectrum. The PACS imaging allows us to 
measure the nebular radius at the wavelengths of the bulk of dust 
emission. To model the dust shell we only considered the spec¬ 
trum at 4 > 20 pm. The spectrum at 4 < 12 pm comes from the 
central star that Voors et al. (12000b fitted with a spherical non- 
LTE model atmosphere. Eor the spectrum between 4 ~ 14 and 
20 pm they argued that it too likely comes from the central star 
and not from some extended source. 

The two-dimensional radiative transfer code 2-Dust (Ueta 
and Meixner 12003b was used to model and interpret the dust 
emission spectrum and the far-infrared images. This is a pub¬ 
licly available versatile code that can be supplied with various 
grain size distributions and optical properties as well as complex 
axisymmetric density distributions. 

It should be mentioned here that since the PACS spectral field 
of view is smaller than the nebular size, the PACS spectrum was 
not taken into consideration for the dust model as this spectrum 
is indeed fainter than the PACS photometric points that contain 
the total flux of the nebula. 

To model the PACS ring of dust with the code 2-Dust, it 
is necessary to consider the morphology of the nebula revealed 
through the infrared PACS and the optical images so as to choose 
the best geometric parameters for the axisymmetric dust density 
distribution model. The 2-Dust code uses a normalized density 
distribution function (Meixner et al. 12002b that is based on a lay¬ 
ered shell model, 

p(R,P)= — 

V-^min/ 

X {l -I- A(1 - cos bY (1) 

where p{R, 9) is the dust mass density at radius R and latitude 
0, T^min is the inner radius of the shell, and Rjw is the superwind 
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Radial distance (") 


Fig. 4. Top left: the 1.5' x 1.5' image of the nebula around AG Car 
observed with PACS at 70 pm. North is up and east to the left. Top right: 
the synthetic image computed with 2-Dust using rjn = 14" and = 
42" and convolved with the PACS PSF. Bottom: Cuts with PA= 160° and 
70° through the central part of the nebula at 70 ^m, observed (black) and 
synthetic (red) ones. 


radius that defines the boundary between the spherical wind and 
the axisymmetric superwind. The first term represents the ra¬ 
dial profile of the spherical wind; the parameters A-F define the 
density profile; the radial factor B can also be a function of the 
latitude through the elongation parameter C; A is the equatorial 
enhancement parameter; the parameter F defines the flatness of 
the shell; and D and E are the symmetry transition parameters 
that describe the abruptness of the geometrical transition in the 
shell. 

It should be specified that we only consider the bright ring- 
nebula for this model. No attempt has been made to model the 
nebular northern extension since it is faint in the infrared and not 
clearly resolved. 

A purely spherical model cannot reproduce the observed 
morphology of the dust shell because there would be too much 
emission at the center of the ring with respect to the observa¬ 
tions. A sphere with equatorial enhancement may be considered 
to reproduce the difference in intensity between the two cuts on 
the infrared image. To appear nearly circular with a clear cen¬ 
tral hole in projection, a spherical shell with equatorial enhance¬ 
ment must be seen at small inclination < 30°, roughly through 
the poles. This scheme is in agreement with the observed global 
morphology. At very low inclination it is difficult to reproduce 
the large difference in intensity between the two cuts in an ax¬ 
isymmetric model like this one. On the other hand, by increasing 
the inclination too much a strongly elliptical ring would be seen, 
which is not observed. In addition to the equatorial enhancement. 



Wavelength {pm) 

Fig. 5. Infrared SED of the nebula around the LBV AG Car from 
data collected at different epochs: the ISO-LWS spectrum and color- 
corrected photometric measurements from IRAS, AKARI and Herschel 
(PACS and SPIRE) space observatories. Results of the 2-Dust model 
fitting are illustrated. Data at T < 20 pm are not considered in the fit. 
The best fit (solid line) is achieved considering two populations of dust 
grains (dashed lines). 


the structure appears clumpy. These inhomogeneities cannot be 
reproduced by the present model. The best approximate axisym¬ 
metric model for the observed ring is a sphere with equatorial 
enhancement (A ~ 5) that is seen at small inclination (~ 30°). 
The values for the other five geometric parameters of Eq. [T]are 
B=3, C=0, D=0, E=0, E=1 to keep the model simple and limit 
the number of free parameters. 

The synthetic image produced by 2-Dust and convolved with 
the PACS PSE is given in the top right panel of Eig.|4] By com¬ 
paring the PACS image to the synthetic one, we determined the 
inner, Hn = 14", and the outer, rout = 42", radii of the dust ring. 
In the bottom panel of Pig. 0] cuts with PA= 160° and 70° through 
the central part of the observed and the simulated nebula are il¬ 
lustrated. Although the basic morphology (radii, thickness and 
axisymmetry of the shell) reproduced by the model agrees with 
the observed one, the intensity of the peaks does not because no 
attempt was made to fit the clumps, as mentioned earlier. 

After constraining the nebular geometry we can proceed with 
the model of the dust SED. Por the stellar parameters we adopted 
the distance D = 6 kpc, the luminosity logL/Lo = 6.1, and the 
temperature - 20000 K (Voors et al. 120^001 Groh et af 12009b . 
The infrared SED of AG Car (Pig. |5]) is too broad to be re¬ 
produced with only one population of dust grains, so we con¬ 
sidered two populations of grains with the same composition 
but different sizes. Voors et al. (12000b showed that the dust in 
this nebula contains large grains, up to 40 pm in radius, and is 
dominated by amorphous silicates with little contribution from 
crystalline species, and more specifically pyroxenes with a 50/50 
Pe to Mg abundance. We therefore adopted a similar dust com¬ 
position. The optical constants of silicates given by Dorschner 
et al. (11995b were used for both populations, extrapolated to a 
constant refraction index in the far-ultraviolet. The size distribu¬ 
tion for the dust grains of Mathis et al. (119771 hereafter MRN) 
was assumed, for each of the two populations: n{a) oc ^ with 
cimm < a < Gimax, where a is the grain radius. The model can be 
adjusted to the data by varying flmax (or amin), which controls the 
20pm / lOO/rm flux density ratio, and the opacity, which controls 
the strength of the emission. 
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Fig. 6. Footprint of the PACS spectral held of view on the image of the 
AG Car nehula at 70 pm. Each number pair is the label of a specihc 
spaxel. North is up and east is to the left. 

The best fit (Fig.|2l was achieved with the use of the follow¬ 
ing populations of dust grains. The first is a population of small 
grains with radii from 0.005 to 1 pm, which is responsible for 
the emission at T < 40 pm. The second is a population of large 
grains with radii from 1 to 50 pm, which is responsible for the 
slope of the SED at d >70 pm. It should be pointed out that 
the large grains are necessary to reproduce the observed infrared 
SED. Several attempts were made to fit the SED using different 
grain sizes. They showed that amax cannot be very different from 
50 pm for the population of large grains. More details are given 
in Appendix A. 

The total mass of dust derived from the modeling is Mdust ~ 
0.21 Mq (0.05 Mq from the small dust grains and 0.16 Mq from 
the large ones), with an uncertainty of ~ 20%. This result is in 
agreement with the total dust mass found by Voors et al. (120001 1. 
The small grains have temperatures from 88 K at rin to 62 K at 
Tout, while the large grains have temperatures from 43 K at to 
29 K at Tout- Similar results are obtained when fitting modified 
blackbody (BB) curves to the SED (see Appendix B for details). 

As described in the introduction, AG Car is a variable star 
with photometric and spectroscopic variations. Although the 
data used for the dust model are obtained at different epochs, 
the visual magnitude of the star is approximately the same at the 
epochs of the IRAS, ISO, AKARI, and Herschel observations. 
In addition, AG Car varies under a roughly constant bolomet- 
ric luminosity (Sect.[T]). By keeping the luminosity constant and 
changing the radius and the temperature of the star within rea¬ 
sonable limits, we do not see significant changes in the model 
results. 

5. Emission iine spectrum 

5.1. Spectrum overview 

Eigure 0 illustrates the footprint of the PACS spectral field of 
view on the image of the nebula at 70 pm. This field of view is 
composed of 25 (5x5) spaxels, each corresponding to a different 
part of the nebula, but it is not large enough to cover the whole 
nebula. 

The integrated spectrum of the nebula over the 25 spaxels is 
shown in Eig.|7] Below 55 pm the shape of the continuum results 


from a spectral response correction in this range that has not yet 
been perfected, while above 190 pm it results from a light leak 
from the second diffraction order of the grating to the first one. 

The following forbidden emission spectral lines are detected 
on the dust continuum; [Oi] AA 63, 146 pm, [Nii] AA 122, 205 
pm, and [C ii] A 158 pm. The absence of higher ionization lines 
indicates that the ionization state of the gas in the nebula around 
AG Car is not as high as in the case of the nebulae around 
other LBVs, for example WRAY 15-751 (Vamvatira-Nakou et 
al. l2013l l. This implies that the gas temperature is lower. 

5.2. Line fiux measurements 

A Gaussian fit was performed on the line profiles to measure 
the emission line intensities in each of the 25 spectra (Fig. |6]) 
using the Image Reduction and Analysis Facility (IRAF, Tody 
I1986ll993ll . The table with these measurements is given in Ap¬ 
pendix C. We note that not all the lines are detected in the outer 
spaxels and that all fluxes reach their highest values at the spaxel 
(3,3). This spaxel corresponds to the southwestern part of the 
nebula, which is the brightest part (see Sect.[3]l. 

Maps of the line intensities were created for each of the five 
detected lines in an effort to investigate whether there are dif¬ 
ferences in the gas properties for distinct parts of the nebula. 
There are only 25 spaxels, and the coverage of the nebula is not 
complete, so we cannot really see the full nebula in these “spec¬ 
troscopic images”. The only wavelength at which we barely see 
the nebular ring is that of 122 pm. Furthermore, maps of line 
intensity ratios of every detected line to the [Nii] 122 pm line 
were created. The differences between distinct nebular regions 
are not significant or convincing. Since these maps are diffi¬ 
cult to interpret given the large uncertainties on the fluxes and 
the wavelength-dependent PSF, they were not considered in the 
present study. 

To measure the total flux of the nebula in each line we need 
to use the integrated spectrum over the 25 spaxels, but as al¬ 
ready mentioned above the beam size of the PACS spectrometer 
is smaller than the size of the nebula. Consequently, the fluxes 
measured using the sum of the 25 spaxels do not correspond 
to the real nebular fluxes. For this reason, the PACS spectrum 
was corrected using the three PACS photometric observations. 
A modified BB was fitted to these data and another one to the 
continuum of the PACS spectrum, for wavelengths smaller than 
190 pm. The ratio of the two curves gives the correction factor. 
In other words, the spectrum was scaled to the photometry. This 
factor linearly depends on the wavelength and goes from 1.11 
at 55 pm to 2.12 at 185 pm. This correction assumes a constant 
line-to-continuum ratio and is analogous to the point source cor¬ 
rection applied in the pipeline to correct the effect of flux lost for 
a point source. 

The [N It] 205 pm line has a problematic calibration in PACS. 
Consequently, its flux needs to be corrected before being used in 
the following analysis. More precisely, the flux in [N ii] 205 pm 
is incorrect for all PACS measurements owing to a light leak, su¬ 
perimposing a lot of flux from 102.5 pm at that wavelength. The 
relative spectral response function (RSRF) used to reduce the 
data suffers from the same light leak. Consequently, when this 
RSRF is applied during the data reduction, the signal at wave¬ 
lengths > 190um is divided by a number that is too high. The 
continuum at these wavelengths is irremediably lost, but pro¬ 
vided one can “scale-back up” with the right number to com¬ 
pensate for the exaggerated RSRF, one can recover the line-flux. 

Using instrument test data obtained on the ground with cali¬ 
bration light sources set at different and known temperatures, one 
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can invert the problem and reconstruct the “clean” RSRF, i.e., in 
the absence of light leak. This suffers from some defects and a 
large uncertainty due to the propagation of errors and to the very 
low response of the instrument at these wavelengths. A correc¬ 
tion factor could nevertheless be derived from it, and confirmed 
within a certain margin by comparison of the line fluxes obtained 
for a few sources by both PACS and SPIRE at that wavelength. 
We finally found that the measured [N ii] 205 pm flux should be 
multiplied by a correction factor of 4.2. An error of 25 % was 
assumed for the final corrected [N ii] 205 pm flux. 0 

The emission line measurements of the nebula integrated 
over the 25 spaxels, before and after the correction for the miss¬ 
ing flux, are given in Table |2] along with the correction factor 
(c.f.) at each wavelength. It should be mentioned that the flux 
measurements of the three lines present in the ISO-LWS spec¬ 
trum of AG Car ([O i] A 63 pm, [N ii] T 122 pm and [C ii] 4 158 
pm) do agree with the corrected values from the PACS spectrum 
within the errors, showing that the correction for the missing flux 
is essentially correct. 


Table 2. Line fluxes of the nebula around AG Car. 


Ion 

A 

(pm) 

F (25 spaxels) 
(10“'5 Wm“2) 

c. f. 

F (corrected) 
(10“'5 Wm“2) 

[Oi] 

63 

7.7 

1.17 

9.0+ 1.8 

[Nil] 

122 

23.6 

1.64 

38.7 + 7.7 

[Oi] 

146 

0.6 

1.83 

1.1 +0.3 

[Cii] 

158 

4.1 

1.92 

7.9 + 1.6 

[Nil] 

205 

1.1 

2.26“ x4.2* 

10.3 + 2.6 


Notes. Missing flux correction PACS/SPIRE cross-calibration fac¬ 
tor 


5.3. Photoionization region characteristics 

The detected emission lines [Nii] 122, 205 pm are associated 
with the Hit region of the nebula around AG Car. The other 
three detected emission lines originate from a region of transi¬ 
tion between ionized and neutral hydrogen and may indicate the 
presence of a photodissociation region (PDR). Their analysis is 
given in the next subsection. 

5.3.1. Hallux 

The Ha-i-[N ii] flux from the nebula was estimated by integrating 
the surface brightness over the whole nebula. Contamination by 
field stars and the background was corrected and the emission 
from the occulted central part extrapolated. The continuum flux 
from the reflection nebula was measured in the adjacent filter, 
accounting for the difference in filter transmissions. However 
since AG Car is a strong emission-line star, the reflected stel¬ 
lar Ha must also be subtracted. Considering the Ha equivalent 
widths measured by Schulte-Ladbeck et al. (I1994I I and Stahl et 
al. (I20011 I for AG Car in 1993-1994 (i.e., accounting for ~ 1.5 
years of time-delay) we estimate the final contamination due to 

^ This part of the infrared spectrum of AG Car has also been observed 
with SPIRE as part of the MESS program. Unfortunately, these data 
cannot be used so as to have a more precise flux for the line [N ii] 205 
pm because the whole ring nebula is outside of the detector coverage 
owing to the geometry of the detector array and because the observing 
mode was a single pointing and not a raster map. Consequently, any 
attempt to recover the nebular flux has huge uncertainty and we decided 
not to include the SPIRE spectroscopic data in our study. 


the reflection nebula to be 20%. The contribution of the strong 
[N It] lines was then subtracted using the [N ii] /Ha ratios from 
available spectroscopic data and the transmission curve of the 
HaH-[Nii] filter. The conversion to absolute flux was done with 
the help of the three spectrophotometric standard stars and three 
planetary nebulae observed in the same filter; the conversion fac¬ 
tors derived from these six objects show excellent internal agree¬ 
ment. 

We measured f’o(Ha) = 1.1 x 10“'° ergs cm“^ s“' uncor¬ 
rected for reddening. The uncertainty amounts to ~ 20%. Adopt¬ 
ing E(B-V) = 0.59 + 0.03 (de Ereitas Pacheco et al. 119921) . we 
derived Eo(Ha) = 4.2 + 0.9 x 10“'° ergscm“^ s ' for the AG Car 
nebula. This flux is higher by a factor of 2 than the fluxes mea¬ 
sured by Stahl (I1987I I. Nota et al. (119921 ) and de Ereitas Pacheco 
et al. (119921 ) in 1986, 1989, and 1991, but is compatible with 
the Hy6 flux measured by Perek (1197H in 1969 (i.e., f’o(Ha) ^ 
3 X 10“'° ergs cm“2 s ' with Ha/H;S = 6 and E(B-V) = 0.59). 
The flux density from the reflection nebula is - 3.9 x 10“'^ 
ergs cm“^ s“' A“' at 6650 A (the central wavelength of the con¬ 
tinuum Alter). The high value of Eo(Ha) we And is in agreement 
with the radio flux, also observed in 1994-1995 (Duncan and 
White I20021 ). and E(B-V) = 0.59 


5.3.2. Electron density 

Smith et al. (119971 ) found a non-constant nebular electron den¬ 
sity, rig, that varies from 600 to 1050 cm“^ using the optical [S ii] 
6731/6717 ratio as an electron density diagnostic. Their result is 
in agreement with those of Mitra and Dufour ( 119901 ) and Nota et 
al. (119921 ) who used the same ratio as a diagnostic. 

In the infrared waveband, the [Nii] 122/205 pm ratio is a 
diagnostic for the electron density of the nebula at low density, 
1 cm“^ < He < 10^ cm“^ (Rubin et al. 119941) . Considering the 
values of Table this ratio is equal to 3.8 +1.2 for the whole 
nebula. The package nebular of the IRAF/STSDAS environment 
(Shaw & Dufour [T9951 ) was used for the calculation of the elec¬ 
tron density. An electron temperature, T^, constant throughout 
the nebula and equal to 6350 ± 400 K was used for all of the fol¬ 
lowing calculations. This is the average temperature calculated 
by Smith et al. (119971 ). The electron density is then found to be 
160 + 90 cm“^. 

The calculated electron density based on the infrared data is 
much lower than the density based on the optical data. This dis¬ 
crepancy is usual and has also been observed in planetary nebu¬ 
lae (Liu et al. 120011 Tsamis et al. 12003b . When the density of a 
nebula is spatially inhomogeneous, different line ratios used as 
density diagnostics lead to different values of the density. This is 
related to the difference in the critical density between the lines 
taken into consideration for the density calculation (Rubin|T989] 
Liu et al. 12001b . The lines [N ii] 122, 205 pm have lower critical 
densities than the lines [S ii] 6731, 6717 A, which means that the 
calculated density using the first pair of lines is smaller than the 
density using the second pair (Rubin fr989b . 

In the following calculations we will use our estimate of the 
electron density based on infrared data because the electron den¬ 
sity is best determined when it is similar to the critical density of 
the lines whose ratio is used as a diagnostic (Rubin et al. 11994b . 
Otherwise, any attempt to calculate ionic abundances will give 
incorrect results (Rubin fr9891 Liu et al. 1200lb . 
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Wavelength (jtim) 

Fig. 7. PACS spectrum of the nebula around AG Car, integrated over the 25 spaxels. Indicated are the lines [Oi], [Nii], and [Cii]. Below 55 jim the 
shape of the continuum results from a spectral response function correction that has not yet been perfected. Above 190 yum the shape results from 
a light leak from the second diffraction order of the grating in the first one. The different observing bands are indicated with different colors. We 
note that the spectral resolution depends on the waveband. 


5.3.3. Ionizing flux 


where Qq (in photons s ') using the Ho' flux is given by 


To calculate the radius of the Stromgren sphere, Rs, and the rate 
of emission of hydrogen-ionizing photons, Qo, a steady nonvari¬ 
able star must be considered. Such an analysis can be done in the 
case of a variable star like AG Car if the recombination time is 
longer than the variability timescale of the ionizing star. The re¬ 
combination time is given by Tree - l/ueO'B yr (Draine 1201 lb . 
where is the recombination coefficient. Using the adopted 
value for the electron density and the assumed electron temper¬ 
ature, the recombination time was estimated to be about 520 yr. 
It is much longer than the timescale of the variability (5-10 yr) 
exhibited by the central star of the nebula, and this conclusion 
still holds if the higher electron density derived in the optical 
is considered. Therefore, an average nonvariable star is a valid 
approximation in our case. 

The values of Qq and Rs can be determined first by using the 
estimated Her flux and second by using the radio flux density, S v 
= 268.7 mJy at 6 cm (4.9 GHz) that was measured by Duncan 
and White (120021 1. adopting a typical error of 0.5 mJy. At 4.9 
GHz the nebula is optically thin and it is assumed to be spherical 
with a uniform density. 

The Rs in pc is given by (Vamvatira-Nakou et al. 12013b 


D Q \ 7 .( 0 . 272 + 0 . 0071 nr 4 ) / Qo 

= liT^J 


( 2 ) 


eo(H«) = 8.59 X lO^^f; (3) 

when using the radio flux it is given by 

eO(radio) = 8.72 X 1043rG0.466-0.02081n7-4) ' ^ ( 4 ) 

where Xe - nejiip is the ratio of the electron to the proton density, 
e is the Ailing factor, T 4 , = TejilO'^ K), v is the radio frequency 
in GHz, D is the distance of the nebula in kpc, Fo{Ha) is the Ha 
flux in ergs cm“^ s ', and is the radio flux in mJy. 

Assuming Xe - 1 because the star is not hot enough to sig¬ 
nificantly ionize He and T 4 - 0.635, we found that the rate of 
emission of hydrogen-ionizing photons is 2 o(h«) = (1-2 ± 0.5) x 
10"^* photons s“' and 2o(radio) = (1-0 ± 0.4) x 10"^* photons s“'. 
There is a good agreement between these two results within the 
uncertainties implying that the value of E(B-V) adopted for 
the calculations is essentially correct. The mean value Qq = 
(1.1 + 0.3) X 10"^* photons s * corresponds to an early-B star 
with Teff ~ 26000 K (Panagia ll973b . which can be considered as 
the average spectral type of the star. 

We also derived Rs = l.l + 0.4 pc assuming e - 1, i.e., that 
the ionized gas Alls the whole volume of the nebula. The fact that 
the nebula is a shell and not a sphere, with inner radius of about 
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Rin — 11" = 0.3 pc in Ha, does not change this result because 
in that case the new Stromgren radius is 
= 1.1 pc. The Stromgren radius is the radius of an ionization 
bounded nebula by definition. In Sect. [2 it was observed that 
the faint part of the nebula in Ha extends up to 0.7 pc from the 
central star. Moreover, the northern faint extension discussed in 
that section extends up to 1 pc. The comparison of these numbers 
with the estimated value of the Stromgren radius, considering 
the uncertainties, leads to the conclusion that the Ha nebula may 
be ionization bounded. The presence of PDR signatures in the 
spectrum supports this conclusion. This value of the Stromgren 
radius is only an average value which can vary locally depending 
on the density inhomogeneities. In particular, according to the 
adopted morphological model, the electron density of the shell 
could be higher along the equator and smaller along the poles 
so that the ionizing radiation can reach the faint extensions or 
bipolar lobes. 


(7) 

In Ha the nebula around AG Car is a shell with inner radius 
Oin = 11 "and an average outer radius 0out = 20". In the radio 
the nebula has approximately the same radii (Duncan and White 
120021) . In this case the mass of the ionized shell nebula is given 
by 

(8) 

The mass of the ionized shell nebula is thus = 6.9+2.8 M© 

and = 6.4 + 2.5 M©, with an average value of = 

6.6 ± 1.9 M©, assuming e - 1. The assumption that the ionization 
of He is negligible (j+ = 0) was made because the central star 
has a temperature lower than 30000 K. This result is slightly 
higher, considering the uncertainties, than the mass of 4.2 M© 
estimated by Nota et al. (119921 [1995]) . 


5.3.4. Abundance ratio N/H 


Given the detected emission lines in the spectrum and the lack 
of [Nm] 57 pm and [Out] 88 pm, only an estimate of the N/H 
abundance number ratio can be made based on the observed Ha 
6562.8 A, [N ii] 122 pm and 205 pm lines and considering that 


N _ <N+) 

H “ ■ 


(5) 


The flux ratios F/Fq(H/T) were calculated for the two infrared 
lines of [N ii] with the observed values of F from Table|2] Using 
the dereddened Ha flux, a case-B recombination with - 6350 
K was assumed to calculate the Hji flux, adopting the effec¬ 
tive recombination coefficient equations of Draine (1201 11) . The 
ionic abundances N^/H^ were then derived using the package 
nebular. The N/H abundance number ratio was calculated to be 
(2.6 ± 1.2) X 10 which is equivalent to a logarithmic value of 
12 + log(N/H) = 8.41 ± 0.20. Considering the errors, this value 
is entirely compatible with that of Smith et al. (119971) . which is 
8.27 + 0.05. It is significantly higher than the solar value (7.83, 
Grevesse et al. 120101) . 


5.3.5. Mass of the ionized gas 


An estimate of the ionized gas mass can be made from the Ha 
and the radio emissions. The equations that are analytically de¬ 
rived in Vamvatira-Nakou et al. (120131) are used for this calcula¬ 
tion. 

For a spherical nebula the ionized mass in solar masses, tak¬ 
ing into account the Ha emission, is given by 


m: 


iTh™ " ^7.9 '^l'+ . 7-^Q-471+QQ151n7'4)^l/2^3/2^5/2pP2^^^', 




(6) 

where 9 is the angular radius of the nebula {R - 6D) in arcsec 
and riii* = «p, ring*, Une** are the number densities of the 
ionized hydrogen, ionized helium, and doubly ionized helium, 
respectively. Assuming nHe++ = 0 and denoting y+ = nHe+/«H+, 
we have Xg = «e/np - 1 + nne*lnn+ = 1 + y+ and p+ ^ 1 + 
4«He+/nH+ = 1 +4y+. 

Considering now the radio flux and using the same formal¬ 
ism as above, the mass of a spherical nebula in solar masses is 
given by 


= 5.82x10-^^^^^ 


+ t-0.175 


_ T 

l+y* " 


/ V nO.05 

Ia9/ 


,1/2^3/2^5/251/2^ 
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5.4. Photodissociation region characteristics 

The fine structure emission lines [Oi] 63, 146 pm and [Cii] 158 
pm are among the most important coolants in PDRs (Hollen- 
bach & Tielens [19971 ). Their detection in our spectrum may in¬ 
dicate the presence of a PDR in the nebula. On the other hand, 
a shock, which is the result of the interaction between the fast 
stellar wind and the slow expanding remnant of a previous evo¬ 
lutionary phase, could also photodissociate molecules and result 
in [O i] and [C ii] emission. However, the values of the calculated 
ratios of [Oi] 63 jum/[Oi] 146 pm and [Oi] 63 / 2 m/[Cii] 158 pm 
are in agreement with the PDR models of Kaufman et al. (119991) 
and not with the shock models of Hollenbach and McKee (11989b . 
In particular, the ratio [Oi] 63 /rm/[Cii] 158 pm is a diagnostic 
between PDR and shock as it is < 10 in PDRs (Tielens and Hol- 
lenbach [l985b . Consequently, based on these ratios, we can con¬ 
clude that a PDR and not a shock is present in the nebula around 
the LBV AG Car and that it is responsible for the [O i] and [C ii] 
emission. Photodissociation regions were detected in the neb¬ 
ula that surrounds the LBV HR Car (Umana et al. 12009b and in 
the nebula around the LBV candidate HD 168625 (Umana et al. 
[20T0b . Later on, the infrared study of the LBV WRAY 15-751 
also revealed the presence of a PDR in the nebula that surrounds 
this star (Vamvatira-Nakou et al. 12013b . 

The physical conditions in the PDR can be determined using 
these three infrared lines, but because of the vicinity of the bright 
Carina nebula, we have to check if these lines come entirely from 
the LBV nebula or if there is a significant contribution to the 
measured fluxes from the background. We therefore checked the 
spectra of the background taken at two different positions on the 
sky and found that the lines [Oi] 63 pm, 146 pm come entirely 
from the nebula. However, the flux of the line [Cii] 158 pm is 
contaminated by the [C ii] foreground/background emission. For 
the nebular spectrum discussed and analyzed in this section, the 
background has been subtracted, as mentioned in Sect.|2] Never¬ 
theless, careful examination of the two off-source spectra shows 
that the background is strong and not uniform. The difference be¬ 
tween the spectra of the two off positions induces an uncertainty 
of at least a factor of 2 on the [C ii] 158 pm line flux. Hence, the 
measured [C ii] 158 pm flux is unreliable and the mass of hydro¬ 
gen in the PDR based on the [C ii] flux cannot be estimated. We 
note that the previous conclusion about the presence of a PDR in 
the nebula is still valid when background contamination is taken 
into account. 

The structure of the PDR is described by the density of the 
atomic hydrogen, and the incident FUV radiation held. Go, 
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which can be calculated using the following equation (Tielens 
12005b . where it is expressed in terms of the average interstel¬ 
lar radiation field that corresponds to an unidirectional radiation 
field of 1.6 X 10“^ erg cm“^ s ', 

Go=625|^. (9) 

where L* is the stellar luminosity, x is the fraction of this lu¬ 
minosity above 6 eV, which is ~ 0.7 for an early-B star (Young 
Owl et al. 12002b . and R is the distance from the star. For the 
PDR of the AG Car nebula, the incident FUV radiation field is 
then Go - 3.7 x 10"', considering that L* = lO^ 'L© (Sect. |4]i 
and R = 0.7 pc, which is the radius of the ionized gas region 
surrounded by the PDR. This result can be used to constrain the 
density of the PDR. The diagnostic diagram of the PDR mod¬ 
els of Kaufman et al. (119991 Figs. 4 and 5) give the ratios of the 
fluxes F[0,]63/7^[Cii]158 and F[Oi]i 45 /F[Oi ]63 as a function of the 
density and the incident FUV radiation field. By using only the 
latter ratio and the calculated Go and considering the uncertain¬ 
ties, we can estimate the density of the PDR to be \ognyp ^ 3, 
with a large uncertainty. 

To verify the consistency of the PDR analysis with the re¬ 
sults of the dust nebula analysis, the dust temperature, Tdust, can 
be estimated based on the radiative equilibrium, since the dust 
absorbs and re-emits the FUV radiation in the far-infrared. In 
case of silicates (i.e., j3 - 2) the dust temperature is given by 
(Tielens l2005l l 

( -j \0.06 ,1/6 

K for Tdust < 250 K . (10) 

We obtain a dust temperature of 7'dust=71 K, assuming a typical 
grain size of a = 0.1 because the average cross-section is 
dominated by small grains. This result is in agreement with the 
results of the 2-Dust model (Sect. S. 

6. Discussion 

The parameters of the LBV AG Car given in Table [2 summarize 
the measurements obtained in this work along with results taken 
from previous studies. The stellar parameters of luminosity, ef¬ 
fective temperature, and distance are from Voors et al. (12000b . 
Groh et al. (12009b . Flumphreys et al. (11989b . and this work. The 
parameters for the shell include the radii, the expansion velocity 
(Smith 11991b . the kinematic age, the ionized gas electron den¬ 
sity and the adopted electron temperature, the abundance ratios 
(N/O from Smith et al. (11997b and N/H from our study), and the 
measured masses of dust and gas. 

The Herschel-PACS infrared images of the LBV AG Car re¬ 
veal a dusty shell nebula that surrounds the central star. It is a 
clumpy ring with an inner radius of 0.4 pc and an outer radius of 
1.2 pc. The Ha-i-[N ii] images show a gas shell nebula that coin¬ 
cides with the dust nebula, but seems to be slightly smaller and 
more elliptical. The nebula has bipolar morphology, a common 
feature among LBV nebulae (Weis l20OT1 Lamers et al. 12001b . 

The nebula around AG Car lies in an empty cavity (Fig. (3). If 
associated with the star, the cavity may correspond to a previous 
mass-loss event when the wind of the O-type progenitor formed 
a bubble, as in the case of WR stars (Marston 11996b . A simi¬ 
lar case is the cavity observed around the LBV WRAY 15-751 
(Vamvatira-Nakou et al. 12013b . though the latter is much larger. 
Velocity mapping of the surrounding interstellar gas would be 
needed to confirm this hypothesis and derive constraints on the 
O-star evolutionary phase. 


Table 3. Parameters of the LBV AG Car and its shell nebula. 


Star 

log L/Lo 

6.1 + 0.2 


Teff (K) 

20000 + 3000 


D (kpc) 

6.0+ 1.0 

Shell 

r\n (pc) 

0.4 


'•out (pc) 

1.2 


Vexp (km S'') 

70 


fkin (10"' yr) 

1.7 


«e (cm'^) 

160 ± 90 


7’e(K) 

6350 ± 400 


N/O 

5.7 + 2.2 


12-(-logN/H 

8.41 + 0.20 


-^dust (M^g) 

0.20 + 0.05 


^^ion.gas (^g) 

6.6+ 1.9 


The results of our study point to a shell nebula of ionized 
gas and dust, surrounded by a thin photodissociation region that 
is heated by an average early-B star. The dust mass-loss rate is 
about (1.8 + 0.5) x 10'^ Mq yr“', considering the duration of the 
enhanced mass-loss episode that was estimated from the kine¬ 
matic age of the inner and outer radii of the nebula. Because we 
do not know the total gas mass as we cannot calculate the neu¬ 
tral gas mass, we must assume a gas-to-dust mass ratio in order 
to estimate the total mass-loss rate. A typical value for this ratio 
is 100 and so the gas mass will be ~ 20 Mq. In the study of the 
nebula around the LBV WRAY 15-751 (Vamvatira-Nakou et al. 
12013b . this ratio was calculated to be about 40. If we assume a 
similar value, the gas mass will be about 10 Mq, higher than but 
comparable to the mass of the ionized gas. Adopting the aver¬ 
age value, the gas mass of the nebula around AG Car is about 15 
Mq with an uncertainty of about 30%. The total mass-loss rate is 
then estimated to be (1.4 + 0.5) x 10“^ Mq yr“'. 

It is interesting to compare this mass-loss rate that corre¬ 
sponds to the period during which the ejection took place with 
recent mass-loss rates. Leitherer et al. (11994b found M{H) - 0.6 
xl0“^ to 4.0 xl0“^ Mq yr“' in 1990-1992 when the star lumi¬ 
nosity was rising, showing no significant dependence on the lu¬ 
minosity phase. Groh et al. (12009b studied the fundamental pa¬ 
rameters of AG Car during the last two periods of minimum, 
1985-1990 and 2000-2003, and calculated a mass-loss rate from 
1.5 xl0“^ to 6.0 xl0“^ Mq yr“'. The mass-loss rate during the 
nebula ejection phase thus appears roughly 50 times higher than 
in the present evolutionary phase. 

The N/O ratio of 5.7 + 2.2 calculated by Smith et al. (11997b 
points to the presence of highly processed material because it is 
much higher than the solar abundances. It is the highest value of 
N/O among the known LB Vs, except the case of rj Car (Smith et 
al. 11998b . The 12-(-log N/H abundance of 8.41 + 0.20, calculated 
on the basis of our observations, is enhanced by a factor of 4.3 
with respect to the solar abundance. It is lower than the value for 
the LBV 77 Car and higher than the values reported for all other 
LBVs (Smith e t al.fTWSl l. 

Groh et al. (12009b calculated the surface abundances of sev¬ 
eral chemical elements at the surface of the star. The comparison 
of the nebular abundances with the surface ones shows that the 
N/O abundance ratio of the nebula is much lower than the surface 
value of 39/jg. As the authors mention, this is compatible with 
the idea that the nebulae around massive stars contain material 
that is less processed than the material of the stellar photosphere. 

Smith et al. (11997b . based on a detailed abundances study, 
argued that the AG Car nebula was formed from material ejected 
during a RSG phase. This was also the suggestion of Voors et al. 
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Fig. 8. Evolution of the N/0 surface abundance ratio as a function of 
the mass-loss rate for a 55 Mg star of solar metallicity and for initial 
rotation rates fi/Ocrit from 0 to 0.4, using the models of Ekstrom et al. 
j2012b . The dashed lines correspond to the adopted value of N/0, with 
its errors, and the lower limit for the mass-loss rate. The thicker lines 
emphasize the part of the tracks compatible with the measurements. Eor 
clarity, the tracks are stopped during the He burning phase (data point 
n“ 195 in Ekstrom et al. 12012b . 

(12000b based on their analysis of the dusty nebula, but Lamers et 
al. (12001b . in their study of the chemical composition of LB Vs, 
concluded that the ejection occurred in a blue supergiant (BSG) 
phase as this can better explain the high expansion velocity. 
Moreover, the problem with an ejection during a RSG phase is 
the lack of luminous RSGs in the HR diagram. 

Based on our observations as well as on evolutionary models 
(Ekstrom et al. 12012b . we can constrain the evolutionary path of 
the central star and the epoch at which the nebula was ejected, 
using the abundance ratios, the measured mass-loss rate, and the 
timescale of the ejection as constraints. The only available abun¬ 
dance ratio that can be used is the N/O ratio. The N/H abun¬ 
dance ratio is indeed sensitive to inhomogeneities of the nebula 
(Lamers et al. 12001b . It should also be stressed that the evolu¬ 
tionary models for massive stars are very uncertain at the post- 
main-sequence phases as they do not include any eruptive event, 
which means that the mass-loss rate recipes are poorly known 
(Smith l20ra . 

A constraint on the initial rotational velocity of AG Car can 
be imposed, based on the results of Groh et al. (1201 lb . In their 
study of AG Car during two periods of visual minimum, they 
concluded that the progenitor did not have a high initial rota¬ 
tional velocity, although they measured the current projected ro¬ 
tational velocity to be 220 km s '. Their conclusions were based 
on the comparison with the evolutionary paths of Meynet and 
Maeder (12003b . The luminosity and effective temperature of the 
star were found to be compatible with the evolutionary tracks of 
a nonrotating star with initial mass between 40 and 60 Mq. 

The total mass-loss rate, estimated during the nebular ejec¬ 
tion, is quite high but uncertain. A lower limit of the mass-loss 
rate can be considered, based on the sum of the dust mass and 
the ionized gas mass that are well determined in the nebula ring. 
Considering the errors, this lower limit is logM = -3.4, where 
M is in Mq yr"'. This result along with the nebular N/O abun¬ 
dance ratio, which is assumed to be the surface abundance ratio 
at the time of the ejection, were compared to the computed evo¬ 
lution of the mass-loss rate versus this abundance ratio using the 
models of Ekstrom et al. (12012b for stars of initial masses that 



Fig. 9. Evolutionary path in the HR diagram of a 55 Mq star of solar 
metallicity and for initial rotation rates fi/Hcrit from 0 to 0.4, using the 
models of Ekstrom et al. (12012b . The thicker lines emphasize the part 
of the tracks compatible with the N/O abundance ratio and the mass- 
loss rate. For clarity, the tracks are stopped during the He burning phase 
(data point n“ 195 in Ekstrom et al. 12012b . 


correspond to the high stellar luminosity of AG Car, considering 
four different cases of stellar rotation from no rotation to a rota¬ 
tion rate of Q/Gcrit=0.4. In Eig.|8] the evolution of the mass-loss 
ratio versus the N/O abundance ratio is illustrated for a 55 M© 
star from the models of Ekstrom et al. (12012b . The measured N/O 
value, with its errors, and the measured lower limit of the mass- 
loss rate are plotted with dashed lines. The part of these tracks 
compatible with the measurements is emphasized with thicker 
lines. To identify at which evolutionary phase of the star this 
corresponds, the same parts of the tracks are reported in the HR 
diagram (Eig. lU). 

Our results are compatible with the evolutionary tracks of 
the models of Ekstrom et al. (12012b for a star of 55 M© with 
solar metallicity and medium rotational velocity. In this case, 
the ejection of the nebula occurs in a post-main-sequence short¬ 
lived episode of high mass loss in agreement with the observa¬ 
tions. We note that in such short-lived episodes, the mass-loss 
rate could be higher than computed from the model since the 
model does not account for eruptive events. Eor a star of 57 M© 
the only compatible evolutionary track is the nonrotating one. 
Consequently, we can conclude that the star may have a low ini¬ 
tial rotational velocity as suggested by Groh et al. (1201 lb . Eor a 
mass of 50 M©, the only compatible evolutionary track is the one 
rotating at Q/Qcut = 0.4. Eor a mass of 60 M©, the N/O ratio is 
reached on the main sequence where the mass-loss rate is much 
smaller than our lower limit, such that no track is compatible 
with both the observed N/O ratio and a short-lived (< 2 10"' yr) 
high mass-loss event. 

A star with initial mass between 40 and 60 M© immediately 
evolves to a BSG without passing through the RSG phase. It 
then evolves towards the LBV and the WR phase (Meynet et 
al. 1201 lb . Groh et al. (12014b performed a detailed study on the 
evolutionary stages of a nonrotating star of 60 M© with solar 
metallicity, combining the evolutionary models of Ekstrom et 
al. (12012b with atmospheric models. Before the WR phase, the 
evolutionary tracks of a 55 M© star (Eig. |3) with little rotation 
are very similar to the track of Groh et al. (12014b for a 60 M© 
without rotation in terms of effective temperature and luminos¬ 
ity. Making use of this result points to an ejection of the nebula 
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during the LBV evolutionary phase of AG Car and more pre¬ 
cisely during a cool LBV phase. Compared to the results ob¬ 
tained for WRAY 15-751, a lower luminosity LBV that passed 
through a RSG phase where the ejection of its nebula took place 
(Vamvatira-Nakou et al. I2013II . this indicates that depending on 
their luminosity, LBV nebulae can be ejected at different evolu¬ 
tionary stages. It should be mentioned that de Freitas Pacheco 
et al. (1 19921 1 compared the AG Car nebular properties, based on 
their spectroscopic observations, with the evolutionary models 
available at that time, and concluded that they were consistent 
with the properties of a star of 60 Mq at the beginning of the 
LBV phase. 

The model of the dust nebula in Sect. |4] showed that large 
dust grains are necessary to reproduce the observed infrared 
SED, in agreement with the results of Voors et al. (12000b . 
This was also the case for the dust nebulae around the LBV 
WRAY 15-751 (Vamvatira-Nakou et al. 12013b and the yellow 
hypergiant Hen 3-1379, a possible pre-LBV (Hutsemekers et 
al. 12013b . Large grains (a > 5 lum) have also been detected in 
supernovae (Gall et al. 12014b . In the case of LB Vs, the stellar 
temperature is most often too high for dust formation to take 
place so that dust production can only happen during large erup¬ 
tions, when a pseudo-photosphere with a sufficiently low tem¬ 
perature is formed. As shown by Kochanek (1201 II 12014b . large 
dust grains can be produced during LBV eruptions and other 
transients when the conditions of high mass-loss rate and low 
pseudo-photosphere temperature are encountered. According to 
these models, to produce grains larger than 10 pm, the central 
star should have gone through a great outburst, with a pseudo¬ 
photosphere temperature as low as 4000 K, i.e., much lower than 
during normal eruptions. During this event, the mass-loss rate is 
expected to be as high as 10“^ MQyr '. For AG Car, this would 
require a duration of the event shorter than estimated from the 
shell thickness, which is possible if the shell thickness is mostly 
due to a spread in velocity (Kochanek l201 lb . 

7. Conclusions 

The analysis of Herschel PACS imaging and spectroscopic ob¬ 
servations of the nebula around the LBV AG Car, along with op¬ 
tical imaging data have been presented. The PACS images show 
that the dust nebula appears as a clumpy ring. It coincides with 
the Ha nebula, but extends farther out. 

The determination of the dust parameters of the nebula was 
performed by dust modeling with the help of a two-dimensional 
radiative transfer code. This model points to the presence of both 
a small and a large grain population of pyroxenes with a 50/50 
Fe to Mg abundance. Large grains (a > 10 pm) are needed to 
reproduce the observational data. 

The infrared spectrum of the nebula consists of forbidden 
emission lines over a dust continuum, without the presence of 
any other dust feature. These lines reveal the presence of ionized 
and photodissociation regions that are mixed with the dust. The 
derived gas abundances show a strong N/O and N/H enhance¬ 
ment as well as a O/H depletion, which is expected for massive 
evolved stars enriched with CNO-cycle processed material. 

The evolutionary path of the star and the epoch at which 
the nebula was ejected were constrained using the abundances, 
mass-loss rate and available evolutionary models. The results 
point to a nebular ejection during a cool LBV evolutionary phase 
of a star with initial mass of about 55 Mq and with little rotation. 
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Fig. A.I. Fits of the far-infrared (A > 20pm) SED of the nebula around 
AG Car. The dust composition (optical constants of silicates with a 
50/50 Mg-to-Fe abundance given by Dorschner et al. 11995b and acm 
are kept constant while the value of changes. 


Appendix A: Tests of the dust model 

In Sect. |4]the model of the dust nebula is presented and the ne¬ 
cessity of a population of large dust grains is stressed. Several 
tests have been performed using the 2-Dust code in an effort to 
reproduce the observed data with different populations of dust 
grains from those adopted in Sect. |4] Our main concern was to 
investigate the influence of the dust grain size and composition 
on the model SED. For this reason, we calculated many models 
by changing one parameter and keeping the other ones constant. 
Keeping in mind that two different populations of grains were 
used to reproduce the broad observed SED in Sect. 0] we con¬ 
sider one population of small grains with radii a^m < a < acen 
and one population of large grains with radii acen < a < ctmax- 
For the dust composition, we used the optical constants given by 
Dorschner et al. (11995b for three different abundances of Mg to 
Fe, 0.5/0.5 (the model presented in Sect. [Hi, 0.4/0.6, and 0.8/0.2. 
For the dust grain sizes, the values of 20, 50, and 80 pm were 
considered for amax and the values of 0.1, 0.3, 1, and 3 pm were 
considered for acen, with the value of amin being kept constant 
and equal to 0.005 pm. 

The comparison of these tests shows that large grains are 
necessary to reproduce the data. The influence of the change of 
flcen on the fit of the observed SED is almost negligible. Fur¬ 
thermore, the fit depends little on the dust abundance of Mg to 
Fe. This is illustrated in Figs. IA.11 - IA.3I In Fig. lA.II three dust 
models are illustrated. The dust composition is the same (optical 
constants of silicates with a 50/50 Mg-to-Fe abundance given by 
Dorschner et al. 11995b and the only parameter that changes is 
Umax- We see that better fits to the data are achieved when large 
grains are considered, in particular when Umax = 50 pm. In all 
cases, when adjusting the observed flux at 250 pm, a value for 
amax that is too small (< 20 pm) gives too much flux at 60-100 
pm. In Fig. IA.2l the comparison of three models with the same 
«max and dust composition but with different flcen is illustrated. 
The model SED depends little on a^en- In Fig. IA.3l the compari¬ 
son of three models with the same grain sizes but different dust 
composition (in terms of the abundance Mg to Fe) is illustrated. 
Again, the model SED depends little on this dust abundance ra¬ 
tio, with a slightly better adjustment for Mg/Fe = 50/50. We fi¬ 
nally note that changing the power law index of the grain size 
distribution does not significantly affect these results. 
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Fig. A.2. Fits of the far-infrared (T > 20/jm) SED of the nehula around 
AG Car. The dust composition (optical constants of silicates with a 
50/50 Mg-to-Fe abundance given by Dorschner et al. 11995b and 
are kept constant while the value of changes. 
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Fig. A.3. Fits of the far-infrared (A > 20/jm) SED of the nebula around 
AG Car. - 50 jum and Ocen = 1 yum are kept constant while the 
dust composition changes. The optical constants of silicates given by 
Dorschner et al. ( 11995b are used for different Mg-to-Fe abundances. 


Appendix B: Modified biackbody fit on the SED 

The observed SED (Fig.|5ll can also be reproduced with a simpler 
model that is the sum of two modified BB curves Fy oc ByCT^y. 
Only the photometric points were considered for making the fit 
illustrated in Fig. ED The mass of the dust can then be derived 
using the equation 


Fig.B.l. Same as Fig.[5]but the infrared SED of the nebula around the 
LBV AG Car is fitted by the sum of two modified biackbody curves. 


second population of dust grains involves large grains (up to 50 
jum) so that this hypothesis is no longer valid. The value of /Cy 
is then taken from the 2-Dust results where it is explicitly com¬ 
puted for the adopted grain population. At 250 /rm, fCy = 13.6 
cm^g which is significantly higher than expected under the 
small-grain approximation. With such large grains the frequency 
dependence of /Cy also differs from the law, being closer to v' 
below 100 yum. A fit with such a composite modified BB gives 
Tdust,i = 76 K and 7’dust,2 = 31 K (which is only slightly different 
from the values obtained with only). Using Edust,! = 76-78 K 
and 7'dust,2 = 31-33 K in Eg. IB.II we then derive Mdust.i = 0.05- 
0.04 Mq and Mdust ,2 = 0.22-0.19 Mq for the small and the large 
grains respectively, in very good agreement with the 2-Dust re¬ 
sults. 


Appendix C: Emission iine fluxes for each spaxei 

The results of the emission line flux measurements for each 
spaxei are given in Table ED The first column contains the de¬ 
tected ions along with the spectral band in which the correspond¬ 
ing line was measured. The following columns contain the line 
fluxes, expressed in W/m^, along with their errors. The spaxei 
numbers (Fig. |6]l are mentioned in every cell of the table. The 
quoted uncertainties are the sum of the line-fitting uncertainty 
plus the uncertainty due to the position of the continuum. 


47dust 


Fy 

By{Ti) Ky ’ 


(B.l) 


where Ky is the mass absorption coefficient, i.e. the absorption 
cross section per unit mass. By the Planck function and D the 
distance to the nebula (Hildebrand 1 19831) . For this calculation, 
the two populations of dust grains are considered independently. 
The fluxes measured at 25 /tm and 250pm are used, i.e., at those 
wavelengths where the contribution of each population domi¬ 
nates (Fig. EU). At 25 pm, Ky - 483 cm^g ' for the silicates 
of Dorschner et al. (119951) . For grains of radii smaller than the 
wavelengths at which dust radiates, Ky is roughly independent 
of the radius and behaves as v® in the far-infrared. However the 
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Table C.l. Line fluxes in each spaxel. A dash indicates a poor S/N or a nondetection. The spatial configuration corresponds to the footprint of the 
PACS-spectrometer as displayed in Fig.[^ 


Ion 

A (band) 

F±AF 

F±AF 

F + AF 

F + AF 

F + AF 


(pm) 

(10-'5 Wm-2) 

(10-'5 Wm-2) 

(10-'^ Wm-2) 

(10-'^ Wm-2) 

(10-'^ W m-2) 



spaxel 0,0 

spaxel 0,1 

spaxel 0,2 

spaxel 0,3 

spaxel 0,4 

[Oi] 

63 (B2A) 

- 

0.40 + 0.05 

0.14 + 0.05 

- 

- 

[Nil] 

122 (RIB) 

0.15 + 0.02 

1.15 + 0.01 

0.85 + 0.02 

0.56 + 0.02 

0.26 + 0.02 

[Oi] 

146 (RIB) 

- 

0.03 + 0.01 

- 

- 

- 

[Oi] 

146 (RIA) 

- 

0.03 + 0.01 

- 

- 

- 

[Cii] 

158 (RIA) 

- 

0.21 +0.01 

0.13 + 0.01 

0.05 + 0.01 

- 

[Nil] 

205 (RIA) 

- 

0.21 + 0.05 

0.11+0.03 

- 

- 



spaxel 1,0 

spaxel 1,1 

spaxel 1,2 

spaxel 1,3 

spaxel 1,4 

[Oi] 

63(B2A) 

0.22 ± 0.05 

0.65 + 0.05 

0.47 + 0.05 

0.29 + 0.05 

- 

[Nil] 

122 (RIB) 

0.64 ± 0.01 

1.86 + 0.02 

1.49 + 0.02 

1.38+0.01 

0.61 + 0.03 

[Oi] 

146 (RIB) 

0.02 ± 0.01 

0.04 + 0.01 

0.04 + 0.01 

0.02 + 0.01 

- 

[Oi] 

146 (RIA) 

0.02 ± 0.01 

0.04 + 0.01 

0.04 + 0.01 

0.03 + 0.01 

- 

[Cii] 

158 (RIA) 

0.09 ± 0.01 

0.25 + 0.01 

0.24 + 0.01 

0.19 + 0.01 

0.06 + 0.01 

[Nil] 

205 (RIA) 

0.10 + 0.02 

0.27 + 0.07 

0.29 + 0.07 

0.24 + 0.06 

0.14 + 0.04 



spaxel 2,0 

spaxel 2,1 

spaxel 2,2 

spaxel 2,3 

spaxel 2.4 

[Oi] 

63(B2A) 

0.41 + 0.05 

0.26 + 0.05 

0.15 + 0.05 

0.65 + 0.05 

0.19 + 0.05 

[Nil] 

122 (RIB) 

1.18 + 0.02 

1.43 + 0.02 

1.01 +0.02 

1.63+0.02 

0.75 + 0.01 

[Oi] 

146 (RIB) 

0.03 + 0.01 

0.02 + 0.01 

0.02 + 0.01 

0.04 + 0.01 

0.02 + 0.01 

[Oi] 

146 (RIA) 

0.03 + 0.01 

0.02 + 0.01 

0 .01+0.01 

0.05 + 0.01 

0.02 + 0.01 

[Cii] 

158 (RIA) 

0.20 + 0.01 

0.24 + 0.01 

0.17 + 0.01 

0.30 + 0.01 

0.14 + 0.01 

[Nil] 

205 (RIA) 

0.17 + 0.04 

0.35 + 0.09 

0.28 + 0.07 

0.33 + 0.08 

0.15 + 0.04 



spaxel 3,0 

spaxel 3,1 

spaxel 3,2 

spaxel 3,3 

spaxel 3,4 

[Oi] 

63(B2A) 

0.25 + 0.05 

0.36 + 0.05 

0.60 + 0.05 

1.35 + 0.06 

- 

[Nil] 

122 (RIB) 

1.04 + 0.02 

1.56 + 0.02 

1.46 + 0.02 

1 .68+0.02 

0.40 + 0.01 

[Oi] 

146 (RIB) 

0.03 + 0.01 

0.02 + 0.01 

0.05 + 0.01 

0.07 + 0.01 

- 

[Oi] 

146 (RIA) 

0.02 + 0.01 

0.03 + 0.01 

0.04 + 0.01 

0.09 + 0.01 

0.02 + 0.01 

[Cii] 

158 (RIA) 

0.16 + 0.01 

0.26 + 0.01 

0.32 + 0.01 

0.43 + 0.01 

0.13 + 0.01 

[Nil] 

205 (RIA) 

0.16 + 0.04 

0.22 + 0.06 

0.23 + 0.06 

0.35 + 0.09 

- 



spaxel 4,0 

spaxel 4,1 

spaxel 4,2 

spaxel 4,3 

spaxel 4,4 

[Oi] 

63(B2A) 

- 

- 

0.58 + 0.05 

0.58 + 0.06 

- 

[Nil] 

122 (RIB) 

0.43 + 0.02 

0.66 + 0.02 

1.06 + 0.02 

0.49 + 0.02 

0.11 +0.02 

[Oi] 

146 (RIB) 

- 

- 

0.04 + 0.01 

0.03 + 0.01 

- 

[Oi] 

146 (RIA) 

- 

- 

0.04 + 0.01 

0.04 + 0.01 

- 

[Cii] 

158 (RIA) 

0.04 + 0.01 

0.09 + 0.01 

0.26 + 0.01 

0.22 + 0.01 

0.04 + 0.01 

[Nil] 

205 (RIA) 

0.09 + 0.02 

- 

0.17 + 0.04 

- 

- 
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